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Abstract 
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tion theory a la Weinberg. The Ml operators have been derived up to N^LO. 
The nuclear matrix elements are evaluated with the use of wave functions ob- 
tained by carrying out variational Monte Carlo calculations for a realistic nu- 
clear Hamiltonian involving high-precision phenomenological potentials like Ar- 
gonne Avl8 and Urbana IX tri-nucleon interactions. We discuss the potential- 
and cutoff-dependence of the results. 
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1 Introduction 



The standard nuclear physics approach (SNPA) based on meson-exchange currents[T], 
high-precision phenomenological potentials and the state-of-the-art techniques for ob- 
taining a few-body nuclear wave functions has achieved tremendous progress in un- 
derstanding light nuclear systems [2]. SNPA is particularly useful when the impulse 
contributions and the well-known long-ranged one-pion-exchanges predominate, as 
was formulated in terms of the so-called chiral filtering [3]. However, if more com- 
plicated meson exchange processes are important, SNPA becomes model-dependent 
because the phenomenological description of the nuclear potentials and transition 
operators does not allow the unique description of their short-range behavior. 

The merit of effective field theory (EFT) [1] is that it offers a systematic way 
to control short-range physics, which in principle allows more accurate theoretical 
predictions than achievable in the conventional approach [5]. An example is the 
hep process, ^He + p — >■ ^He + Ue + e"^ , which is potentially important for solar 
neutrino physics. A reliable estimation of the S'-factor for this reaction has been a 
long-standing challenge in nuclear physics [6]; its theoretical estimates have varied 
by orders of magnitude in the literature. A highly elaborate SNPA calculation of 
this S-factor was carried out by Marcucci et. al. but it was a recent EFT-based 
calculation by Park et al. [S] that provided definitive support to the SNPA results, 
giving in addition a quantitative error estimate. In [8], the authors have developed 
an EFT approach, called EFT* or MEEFT {more effective effective field theory) and 
performed a parameter-free calculation of the hep cross section, which has led to an 
estimate of the S'-factor with ~15% precision; for a detailed review, see Ref. [9]. The 
same method has marked successes also in describing the highly suppressed isoscalar 
amplitude pertinent to the n+p— process [IOl[TT], as well as the weak processes 
like muon capture on the deuteron [12] and v-d scattering [T3l [T3]. For a recent 
review, see Ref. [15]. 

The strategy of MEEFT, as explained in [3 El [IS], is in the spirit of the original 
Weinberg scheme [T7] based on the chiral expansion of irreducible terms. In MEEFT 
the relevant current operators are derived systematically by applying HB^PT to a 
specified order. Nuclear matrix elements corresponding to these current operators are 
evaluated with the use of reahstic nuclear wave functions obtained by applying an ab- 
initio or quasi-ab-initio few-body calculation method to a realistic phenomenological 
nuclear Hamiltonian that involves high-precision phenomenological potentials. In 
MEEFT, as in the usual EFT, short-range physics is simulated by contact counter- 
terms and the coefficients of these terms, called the low-energy constants (LECs), are 
determined by requiring that a selected set of experimental data is reproduced. This 
renormalization procedure is expected to remove, to a large extent, model-dependence 
that might creep in through the phenomenological parametrization of short-range 
physics in the adopted nuclear interaction. We note that the strategy employed 
here is closely related to the one used to construct a universal Viow-k in the recently 
developed renormalization-group approach to nuclear interactions |18j . 

The purpose of this letter is to demonstrate that MEEFT provides a practical 
and reliable tool to compute the low energy Ml properties of few-nucleon systems 
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in a largely model-independent way. We calculate here four low-energy electromag- 
netic observables, namely the magnetic moments of ^H, ^He and and the rate 
of radiative capture of a thermal neutron on a proton. We derive all the relevant 
operators up to the next-to-next-to-next leading order (N^LO) of the perturbation 
series, which include the short-range contact terms. We use the experimental values 
of two of the above-mentioned four observables to fix the coefficients of the contact 
terms appearing at N^LO and make prediction for the remaining two observables. 
We test the consistency as well as the model-independence of the obtained results. 
To our knowledge, this is the first EFT calculation of the magnetic moments of the 
A=3 systems. 

Before closing this section, we remark that other EFT methods can also be 
adopted, at least in principle, for studying the Ml properties of few-nucleon systems 
at low energy. One of such possible alternatives consists in deriving not only the tran- 
sition operators but also the nuclear interactions using the same EFT. This method 
has the advantage of being more transparent in discussing the order counting and 
model-independence; for the recent developments in the EFT studies of nuclear po- 
tentials, see Refs. [19], [201 [21] and references therein. Such a study is expected to be 
highly illuminating particularly for the electric transitions because gauge invariance 
(or charge conservation) is automatically guaranteed order by order. Another possi- 
bility is to use the so-called pionless EFT, where even the pions are integrated out 
leaving only the nucleon field as a pertinent degree of freedom. Pionless EFT has been 
successfully applied to some processes involving two- or three- nucleon systems [22] . 

2 Current operators 

The relevant operator for the magnetic moment and radiative np capture at 
threshold is the Ml operator n{q), 



where T^Q°'^{q) is defined in [23], = {u, q) is the momentum carried out by the 

photon (g^ = for the magnetic moment while g'^ 7^ for np capture), and q= \q\. 

We derive the Ml operator in HBxPT, which contains the nucleons and pi- 
ons as pertinent degrees of freedom with all other massive fields integrated out. In 
HBxPT the electromagnetic currents and Ml operator are expanded systematically 
with increasing powers of Q/A^, where Q stands for the typical momentum scale 
of the process and/or the pion mass, and ~ 4:'^f-K ~ m ~ 1 GeV is the chiral 
scale, ~ 92.4 MeV is the pion decay constant, and m is the nucleon mass. We 
remark that, while the nucleon momentum is of order of Q, its energy (~ pf/m) 
is of order of /m, and consequently the four-momentum of the emitted photon q^ 
(with uj = \q\) should also be counted as 0{Q'^/m). In this work we include all the 
contributions up to N'^LO, where N'^LO denotes terms of order of {Q/A^Y compared 
to the leading one-body contribution. It is worth mentioning that there is a different 




(1) 
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power counting scheme where the nucleon mass is regarded as heavier than the chi- 
ral scale, see Refs. [191 120] for details. However, the use of this alternative counting 
scheme would not affect the results to be reported in this article since the difference 
between the two counting schemes would appear only at orders higher than explicitly 
considered here (N^LO). 

The one-body (IB) Ml operator including the relativistic corrections reads 



(^i ( fj'i — Qi 



2m^ 



l-^i Qi — — 
O 2 Pi^i ■ Pi 



+ 



30 



r^ - ( ^ P^ \ w{2fii - Qi) . ' 

""^^y 2^) ^ — '''' ""'^ 



(2) 



where j„(x) = ^^"^t^^'"' jn(a^) = 1 + 0{x'^), Qi = {1 + Tf)/2 is the charge of the i-th 
nucleon, and = (/is + rf/i„)/2 is the magnetic moment in units of the nuclear 
magneton, fi^ = e/(2mp), with Us = l^p + jJ'n — 0.8798 and = Hp — Hn — 4.7059; 



2V V I "i Vt) with the understanding that the derivatives act only on the wave 
functions, and = |rj|. In the above equation the familiar (/ijCTj + QiVi x term is 
of leading order (LO), while all the others terms are N^LO. Since there are no N^LO 
contributions to the H^b-, the neglected terms are of N'^LO and higher orders. 

Corrections to the IB operator are due to the meson-exchange currents (MEC). 
Up to N^LO, only two-body (2B) contributions enter; three-body (3B) currents are 
N^LO or higher order. It is to be emphasized that MECs derived in EFT are mean- 
ingful only up to a certain momentum scale A, where A is the cutoff below which the 
chosen explicit degrees of freedom reside. This cutoff may be realized by introducing 
a Gaussian regulator in performing the Fourier transformation of the MECs from mo- 
mentum space to coordinate space [8]. It is to be noted that the contributions due to 
high momentum exchanges (above the cutoff scale) are not simply ignored but, as we 
will discuss later, they are accounted for by the renormalization of the contact-term 
coefficients. 

We decompose the two-body current into the so/i-one-pion-exchange (Ivr) term, 
vertex corrections to the one-pion exchange (IttC) term, the so-called fixed term 
contribution (Ivr : fixed), the two-pion-exchanges (2tx) term, and the contact-term 
contribution (CT), 



I \ r Itt , / IttC , ln:fixed , 2n , CT 



NLO + N^LO. 



(3) 



i<j 



The leading MEC due to the so/t-one-pion-exchange (Ivr) is NLO and given as 



P'12 



24/2 



T, 



(X) 



MqR) 



T'^Rxr[(Ti- (T2y^^{r) + Si2y2A{r)]UqR) 



(4) 
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(a) (b) (c) (d) (e) 

Figure 1: Tree graphs of NLO and N^LO. One-pion exchange "seagull" (a) and "pion- 
pole" diagram (b) contribute to the fi^'^. Diagrams (c)-(e) contribute to the fi^'"^ and 
^injixed N^LO. The dot represents the vertex corrections coming from NLO or 
N^LO lagrangian. 
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where qa — 1.2695, r = ri — r2, r = \r\, R = {ri + ^2)72, R = \R\, and 
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(27r) 

^lA = -^jfVoA^ ylk = ^i^l^i^yoA and S12 = 3cri ■ fcrs ■ f - cri ■ cts. We have also 

^2, 



defined T^®^ = TqCTq and T^®^ = r| [f f ■ ctq — Ictq] , = Ti r2, ctq = cti 0-2 
where = ±, x . Strictly speaking, contains not only NLO but also admixtures 
of N^LO and higher order terms; if we expand fi^J in powers of q, the g-independent 
part is NLO while the remaining g-dependent part is N^LO (or of even higher orders). 
Furthermore, if q is non-zero, as is the case for the np capture, Fourier transformation 
for the so/i(-one-pion-exchange current becomes rather involved, generating compli- 
cated N^LO contributions, which are denoted by the ellipsis in the above equation and 
will be reported elsewhere |2l]. A numerical evaluation shows that the contribution 
of the N^LO part in the Itt current is negligibly small. 

There are no corrections at N^LO, and all the contributions up to N'^LO are 
included in this work. Among the N^LO contributions, the one-loop vertex correction 
of the one-pion exchange has been investigated in detail in Ref . [TOl [TTj , 



t^lf = -||(c. + ca) 



jiiqR) 



9a 



1 ^ 1 ^ 

g-^S yOA 2^ 



NwzTi ■ T2 [cr+?/oA + (St-^- ■ cr+ - cr+)y^J ]i{qR)- (6) 



16/r 

The values of the LECs, (c^^, ca, -^vyz), should in principle be fixed either by solving 
the underlying theory, QCD, or by fitting to suitable experimental observables. Since 
this has not yet been done, we adopt here the estimates given in Ref. [TOlfTT] based on 
the resonance saturation assumption and the Wess-Zumino action, (c^, ca, Nwz) — 
(0.1021, 0.1667, 0.02395). Further discussion on the use of these estimates will be 
given later in the text. 

The "fixed-term" contributions, nlj'^^^'^'^, represent vertex corrections to the 
5o/t-one-pion-exchange and fixed completely by Lorentz covariance. They might be 
viewed as relativistic corrections to the irNN and A'^^irNN vertices. These correc- 
tions can be obtained conveniently by performing the Foldy-Wouthuysen transforma- 
tion of the relativistic Lagrangian, and the resulting vertex functions read 
• 7i{k) + N{p) -> N{p') vertex : 

+ -^^[4:p ■ k(T ■ p- Ap^o- ■ k-2p - kcr ■ k + 2k^(T ■ p - k'^cr ■ (7) 
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71" (k) + N{p) N{p') + (g) vertex : 

^l^e^-^r^l^ (l-^ [V -Ap.iq-k) + 2k'-2k.q + q']^ 

4m2 ' P ^ P^ ' P + ' C^l ~ ^) ~ ■ k + ik X J. 

S^-)l}^fj + J_[kx {p X a) + a X {p X k) + tq X k]}. (8) 
1 2m 4m"^ J 



} 



2// 



The expressions of the Ml operator corresponding to the fixed terms are quite lengthy 
and will be reported elsewhere [21]. The fixed terms containing the nucleon momen- 
tum operators make the calculation highly involved. 

The two-pion-exchange diagrams shown in Figj2]give rise to 



^-12 



1 



Ls + ( 4"-^ 



— ^(ri X T2yR X r—Lojo{qR) 



f^-^ ) L 



2bQTi^n 



dr 



(9) 



where 



f r|-J^o + 4 [4^1 - + r^SKo + 2K, 



dr 



2K2 + 9l{8K2 + 2Ki + 2Ko) 



(10) 



The loop functions K's are defined in Ref. [H [TO] . 

Finally there are contact-term contributions of the form 



2m, 



:iii 



where (745 = 771^(74 and (74^ = — mp(G'^ + ^E^'^)] (74, G;^ and E^'^ are the coefficients 
of the contact terms introduced in Refs. pHl [TT] . A noteworthy point is that, after 
removing redundant terms, there are only two independent contact-terms relevant 
for the Ml operator up to N^LO. This reduction of the effective number of counter 
terms is due to Fermi-Dirac statistics; a similar reduction has been noticed for the 
Gamow- Teller operator [8], where only one linear combination of LECs needs to be 
retained. 



3 Results 

We now use the Ml operators described above to calculate the magnetic mo- 
ments of '^He, and as well as the rate of radiative capture of a thermal neutron 
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on a proton, np d^y. The total np — > d'j cross section at threshold reads 



np 




(12) 



with 



Mnp = 2mp / d^riplQ{r)tJ,^{q)iJnp{r) 



(13) 



where Vn = 2,200 m/s denotes the neutron velocity in the lab frame, ipd,o{r) the 
deuteron wave function with the spin component = 0. The ipnp{i^) denotes the 
spin-singlet np scattering wave function obtained by applying the normalization con- 
dition lim^^oo Mnp(r) = r — as, where is the spin-singlet np scattering length. The 
experimental value of the radiative np capture cross section is 0"^^^ = 332.6 ± 0.7 
mb [25], which corresponds to M^^p = 410.2 ± 0.4 fm^^l 

A realistic nuclear Hamiltonian is constructed using a high-precision phenomeno- 
logical potential. We consider here the Avl8 bi-nucleon potential [26] with/without 
the Urbana IX (U9) tri-nucleon force (3NF) [27], and we also consider the Argonne 
Avl4 bi-nucleon potential with/without the Urbana VIII tri-nucleon force (US). We 
remark that, from the HBChPT point of view, 3NF is N^LO compared to the leading 
nucleon-nucleon interactions, and hence should be included in our N^LO calculation. 

The variational Monte Carlo (VMC) technique [271 [28] is used to solve the 3- 
body Schroedinger equation in order to obtain wave functions of and ^He. MEEFT 
is based on the assumption that the adopted nuclear wave functions are exact to 
the order under consideration. This presupposes (i) the validity of using the high- 
precision phenomenological potential instead of an EFT-based potential, and (ii) the 
sufficient accuracy of the method used for solving the A-body Schroedinger equation 
for a given phenomenological A-body Hamiltonian. The first point is a fundamental 
issue that awaits further detailed studies. As for the second point, we remark that 
the long-range behavior of the wave functions can significantly affect the calculated 
nuclear transition matrix elements, and hence one should avoid using schematic wave 
functions. Unfortunately, it is difficult to fully quantify the uncertainty related to 
the wave functions used in this work. One of indirect measures for the accuracy of 
the adopted potentials and the VMC technique is the binding energies of and 
^He. We give in Table [1] the results obtained in VMC for the AvlS and Avl8-hU9 
potentials; also shown in the table are the results obtained by solving the Faddeev 
equations [291 EOl EI]. Note that the use of 3NF is imperative in order to reproduce 
the binding energies of and ^He. Later in the text we shall give comparison of 
the results for the Ml observables calculated with the Avl8(-|-U9) potential and the 
Avl4(-|-U8) potential, and this comparison provides another indirect measure for the 
stability of our results against changes in the input potential. 

The cutoff A has a physical meaning, and its choice is not arbitrary [32]. Thus A 
should be smaller than the masses of the vector mesons that have been integrated out. 
Meanwhile, since the pion is an explicit degree of freedom in our scheme, A should 
be much larger than the pion mass in order to ensure that all pertinent low-energy 
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Table 1: The binding energies (in MeV) of and ^He calculated in the VMC method 
for the Avl8 and Avl8+U9 potentials. For comparison, we also give in the square 
brackets the results obtained by solving the Faddeev equations. 





Avl8 


Avl8+U9 


Exp 


BE(3H) 


7.35(1) [7.61] 


8.24(1) [8.47] 


8.48 


BE(3He) 


6.59(1) [6.91] 


7.48(1) [7.74] 


7.72 



Table 2: Magnetic moments of and '^He in units of the nuclear magneton calcu- 
lated for the Avl8+U9 interactions. Also listed are the values of g^s and g^^ which, 
for a given value of A, reproduce the experimental values of /i(^H) and Mnp- The 
parenthesized numbers represent the Monte Carlo statistical errors, while the en- 
tries in the square brackets correspond to the calculations in which the contact term 
contributions are ignored. 



A [MeV] 


^(^H) 




94s 




500 


3.034(13) [2.883] 


-2.196(13)[- 


-2.074] 


0.5786 


2.8790 


600 


3.035(13) [2.944] 


-2.198(13)[- 


-2.120] 


0.2995 


1.9567 


700 


3.036(13) [2.988] 


-2.199(13)[- 


-2.150] 


-0.0202 


1.2954 


800 


3.037(13) [3.019] 


-2.200(13)[- 


-2.168] 


-0.3965 


0.7882 
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Table 3: The magnetic moments of ^He, '^H, (in units of /itv) and the np d'j 
matrix elements (in units of fm^''^) calculated with the Avl8+U9 potential for A = 600 
MeV. 





u(2H) 




/i(3H) 


u(3He) 


LO: IB 


0.8469 


393.1 


2.585 


-1.774 


NLO: 


Itt 


0.0000 


8.7 


0.205 


-0.205 


N^LO 


IB l/m2 


-0.0069 


-0.1 


-0.018 


-0.007 


N^LO 


IttC 


0.0077 


4.4 


0.133 


-0.116 


N^LO 


fixed term 


0.0044 


-0.2 


-0.003 


0.014 


N^LO 


2tt 


0.0000 


1.5 


0.043 


-0.043 


N^LO 


g4^s term 


0.0053 


0.0 


0.007 


0.007 


N^LO 


g^^ term 


0.0000 


2.9 


0.085 


-0.085 


N^LO total 


0.0174 


8.6 


0.265 


-0.223 


total 


0.8574 


410.2 


3.035 


-2.198 


experiment 


0.8574 


410.2(4) 


2.979 


-2.128 



contributions are properly included. In the present work, we consider A = 500, 600, 
700 and 800 MeV as representative values. 

For a given value of A, we adjust (74^ and Qa^v to reproduce the experimental 
values of /i(^H) = 0.8574/^7v and (or equivalently M^^). We remark that the 
fitted values of (74^ and g^^ depend on a particular choice of the potential model used 
to calculate the wave functions as well as the cutoff parameter. In Table [2], we show 
our predictions for the and ^He magnetic moments obtained with the Avl8+U9 
potential for the four representative values of A; also shown are the values of g^s and 
(74„ optimized for each value of A. To highlight the roles of the contact terms, we 
include in the table the results corresponding to cases for which the contact terms are 
artificially dropped; see the entries in the square brackets. Table [2] indicates that the 
results of the full calculation are almost independent of A, whereas those obtained 
without the contact-term contributions show pronounced cutoff dependence. 

Table [3] shows the contributions of the individual Ml operators to the magnetic 
moments of ^He, and as well as the matrix element for np d'j, calculated 
with the Avl8-MJ9 potential for A = 600 MeV. The table indicates that the NLO Ivr 
contribution is rather small, comparable to the total N^LO contribution. Substantial 
cancellation among the various terms contributing to fi^'^ is responsible for this fea- 
ture. Thus the /j,^'^ contribution to M„p, 8.7 in units of fm^''^, can be decomposed as 
8.7 = (19.4 — 14.4 -|- 4.0) — 0.3, where the first three terms come from the first line of 
eq.(jl]), and the fourth term from the second line. 

Finally, in Table H] we compare the results for the different choices of the nu- 
clear potentials; the table also gives the results of the SNPA calculations taken from 
Refs. [33| [M] , as well as the experimental data [35] . The numbers in the parentheses 
represent statistical errors of our VMC calculation. We remark that the uncertainty 
in the experimental value of (T„p, which affects the determination of g^s and(y'4„, causes 
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Table 4: The A = 3 magnetic moments (in units of /iat) calculated for A = 600 MeV. 
The last column gives the values of g^s and g^y fitted to reproduce the experimental 
values of /i(^H) and M„p. The bottom three rows show the SNPA results and the 
experimental data. 





^(3H)+^(3He) 


/i(3H) - ^(^He) 


94:S gAv 


This work (Avl8+U9) 


0.838(0) 


5.233(25) 


0.300(6) 1.96(29) 


(Avl8) 


0.838(0) 


5.242(26) 


0.300(6) 1.96(29) 


(Avl4+U8) 


0.838(0) 


5.266(33) 


0.391(7) 2.25(31) 


(Avl4) 


0.844(0) 


5.204(30) 


0.391(7) 2.25(31) 


SNPA I [33] 


0.828 


5.078 




SNPA II [31] 


0.884 


5.114 




Experiment |35j 


0.851 


5.107 





~ 0.5 % uncertainty in yu(^H)-yu(^He). Our results are found to be almost independent 
of the interaction model and the cutoff values used, satisfying the general tenet of 
EFTs that the model- and cutoff-dependence should be of higher order than the order 
to which calculation is performed, i.e., N^LO for the present case. This potential- 
independence is also consistent with the notion of the Viow-k [3S] , which dictates that 
all the high-precision phenomenological potentials become universal if the cutoff is 
lowered to 2 fm~^. Physically, this reflects the fact that the differences among various 
high-precision potentials lie only in the high-energy region that can be renormalized 
away leaving little effect for low-energy dynamics, provided that the renormalization 
procedure is correctly done. 

Comparison of the results with and without the tri-nucleon forces (3NF) in 
Table H] indicates that the role of 3NF in the magnetic moments is small; the largest 
effect seen for the Avl4-|-U8 case is of the level of 1%. 

Our no-parameter MEEFT calculation gives the values of /i(^H) + /i(^He) and 
//(^H) — /x(^He) which agree with the experimental values at the 2-3 % level. This 
is slightly worse than the agreement reported in the SNPA calculations. We will 
come back to this point in the next section. We emphasize, however, that as far as 
the structure of the Ml operators is concerned, the present work gives a complete 
expression, a feature that distinguishes MEEFT from SNPA. 

4 Discussion 

We have reported here the calculation of the magnetic moments of and ^He 
based on HBChPT. All the Ml operators up to N^LO have been explicitly derived. At 
N^LO, two unknown parameters, (74^ and g^v, enter as the coefficients of contact terms. 
Following the MEEFT strategy, we have fixed them by imposing the renormalization 
condition that the experimental values of the deuteron magnetic moment and thermal 
neutron capture rate on proton be reproduced at N^LO. 
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Concerning the 2 ~ 3 % discrepancy between the values of /i(^H) + /i('^He) and 
/i(^H) — /i('^He) obtained in our MEEFT calculation and the corresponding experi- 
mental values, there are a few points to be discussed. The first is the level of accuracy 
of the A = 3 wave functions obtained in the VMC method. The indication that the 
VMC wave functions deviate slightly from the accurate wave functions is already visi- 
ble in Table [1], which shows that the binding energies of the A = 3 systems calculated 
in VMC do not quite agree with those obtained in the Faddeev calculation [29 | [3 m [3T]. 
An attempt is in progress [21] to improve the present work with the use of ab initio 
wave functions obtainable in the Faddeev method and to extend our formalism to the 
^He 4- n — > ^He + 7 reaction. 

The second point is higher order contributions. As can be seen in Table 3, the 
difference between theory and experiment is only about 30 % of N^LO contributions, 
for both isovector and isoscalar channel. Given the HB^PT expansion parameter, this 
suggests that N^LO effects of natural size can remedy the discrepancy. Among the 
N^LO contributions, we expect the three-body (3B) currents - which appear first at 
N^LO and hence have not been studied here- to play the most important role for the 
following reasons. We note that a substantial portion of the higher-order two-body 
contributions can be effectively absorbed into the renormalization of the coefficients 
of the contact terms; this is particularly true for short-ranged contributions. To illus- 
trate this aspect, we consider the role of the 27r contribution in our N'^LO calculation 
and compare the results of two "27r-contribution-less" calculations. In the first case, 
the 27r terms are simply dropped without any other accompanying changes, whereas 
in the second case the values of g^s and g^^ are readjusted to reproduce the experi- 
mental values of /i(^H) and (T„p {without the 27r contributions). We expect that the 
change in the net results for the second case, which may be viewed as an effective 2tt 
contribution, is much smaller than the 27r contribution estimated from the first case. 
We have verified this feature numerically. From this example, we may expect that 
the effective N^LO 2B contributions are small. The situation is however completely 
different for the 3B contribution in the Ml operator, Since /Xgg appears in the 

^ = 3 systems but not in the A = 2 systems, its contribution cannot be absorbed into 
the two-nucleon contact terms. The higher order calculation up to N^LO is relegated 
to future work. 

Finally, we discuss the uncertainty related to the values of the LECs, c^, ca and 
N\Yz, which appear in /xjg*", eq. ([6]). The /xjg*" contributions themselves are compa- 
rable in size to the difference between the present calculation and the experimental 
data, as can be seen in Table O However, due to the above-mentioned renormaliza- 
tion procedure, the effective contributions are expected to be small. We have 
checked this by turning off ^^12*" and re-performing the renormalization procedure, 
and found that the effective ^12*^ contributions are negligible. If on the other hand 
we naively treat these LECs as free parameters and adjust them to reproduce the ex- 
perimental values, we are led to very large values of the LECs, which strongly violates 
the naturalness condition. 
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